Emergence of a negative charging energy in a metallic dot capacitively coupled to a 

superconducting island 
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We consider the hybrid setup formed by a metallic dot, capacitively coupled to a superconducting island 
S connected to a bulk superconductor by a Josephson junction. Charge fluctuations in S act as a dynamical 
gate and overscreen the electronic repulsion in the metallic dot, producing an attractive interaction between two 
additional electrons. As the offset charge of the metallic dot is increased, the dot charging curve shows positive 
steps (+2e) followed by negative ones (— e) signaling the occurrence of a negative differential capacitance. A 
proposal for experimental detection is given, and potential applications in nanoelectronics are mentioned. 
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PACS numbers: 73.23.Hk, 74.78.Na 



At low temperatures, the electronic transport through point- 
like metallic nanostructures (quantum dots) is dominated by 
the electronic Coulomb repulsion between additional elec- 
trons. When a small-capacitance dot is weakly coupled to 
a normal metallic reservoir, the average number of charges 
in the dot, n^, increases one by one with the gate voltage 
V 9 n, leading to conductance peaks [1]. This Coulomb block- 
ade phenomenon has recently enabled an individual control of 
charge or spin, for instance in view of quantum information 
protocols (|2|, |3[] . We address here the possibility of inverting 
the sign of the charging energy. Indeed, creating a negative 
effective charging energy in a normal (non-superconducting) 
metallic dot would induce attractive correlations, triggering 
for instance pair tunneling from/to a normal reservoir ||4j] or a 
charge Kondo effect [5, 6], or giving rise to super-Poissonian 
shot noise |@, 01 . Going beyond the single dot case, attrac- 
tive correlations between electrons in spatially separated dots 
could help implementing quantum information protocols in- 
volving two-qubit gates, when the qubits are carried by the 
charge (spin) of the last added electron. Historically, attractive 
interactions in the solid state are known as valence-skipping 
states [8], and negative-U centers 091]. Another possible mech- 
anism for electronic attraction is mediated by optical phonons, 
binding two electrons as a bipolaron in confined geometries 
for strong electron-phonon coupling and lattice polarizability 
ifToh . Due to the low polarizability and small effective carrier 
mass, bipolarons are unlikely to form in a clean GaAs/AlGaAs 
two-dimensional electron gas (2DEG), although they might 
do so in presence of a few donor impurities fiUl . or with 
more polarizable multilayer materials 1131 . Molecular junc- 
tions are also promising for achieving a negative charging en- 
ergy HQ. 

In the present Letter, we propose an alternative mechanism, 
pointing out that the repulsive charging energy in a metallic 
dot (N) connected to a normal reservoir (Fig. 1) can be turned 
into an attractive one when N is capacitively coupled to a su- 
perconducting island (S). The latter is connected to a super- 
conducting reservoir by a Josephson junction (J J) and operates 
in the Cooper pair box regime, e.g. it fluctuates between two 
pair number states 014ill5ll . Here we assume that electron tun- 



neling between S and N is negligible; therefore no proximity 
effect occurs in the N dot. We instead focus on the charging 
properties of the N grain as its gate voltage is varied. The S is- 
land acts as an effective dynamical gate, whose effects turn out 
to be nonlinear. As the main result of this Letter, the Coulomb 
charging energy in N can be overscreened by the neighboring 
pair fluctuations in S, and an effective local attraction appears 
between electrons added into N. As a corollary, certain charge 
states are "skipped" as the N gate voltage is varied. The re- 
sulting charging curve becomes non-monotonous, displaying 
positive steps (+2e) followed by negative ones (— e). A re- 
lated screening effect was proposed by Averin and Bruder for 
controlling the coupling between two superconducting charge 
qubits [16]. Notice that if N were coupled to both drain and 
source reservoirs, our set-up would be similar to a Cooper pair 
box coupled to a single-electron transistor (SET). The latter 
has been studied in great detail as a read-out device for a su- 
perconducting (charge) qubit embodied in the S island lfl7ll . 
In this case, the capacitive coupling between N and S must be 
small, in order to minimize the decoherence due to backaction 
of the normal part of the device onto the superconducting one, 
whereas in our proposal, the coupling is very strong. 
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FIG. 1 : Schematic view of a normal grain (N) coupled with a strong 
capacitive coupling (controlled by 3, 9) to a Cooper pair box com- 
posed of a Josephson junction connecting superconducting reservoir 
2 and island, and gate 10. Electrons tunnel between N and its reser- 
voir (defined by 7, 8). Detection is made by sweeping the gate volt- 
age 4 and measuring the island voltages using quantum point contacts 
for both N (5,6,7) and S (1,11,12). 
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The JJ connecting the S island to the reservoir has a Joseph- 
son energy Ej and capacitance Cj, and a gate imposes an 
offset vg = CgsV g s/e, with C g g < Cj. Similarly, the N 
island is connected to a normal reservoir by a tunnel junc- 
tion, with single-electron tunneling rate T and capacitance 
Cn, and experiences a gate offset vn — CqnVqn /e, with 
C g N <C Cn- Most importantly, the islands N and S are 
coupled by a large capacitance Co > Cn, Cj. We assume 
the superconducting gap in S to be larger than the charg- 
ing energy, such that only even charge number states ng oc- 
cur in S. At low temperatures, quasiparticle tunneling in S 
can be neglected. Defining C^s — Cj + Co + C g s and 
Csn = Cn + Co + C g N, and introducing the parameters 
b = Csat / Css and r — Co/ %/CsatCss, the total charging 
energy of the NS system can be written as 111 811 

E c = E C N[{n N - v N ) 2 + b(ng - v s ) 2 

+ 2rVb(n N - v N )(ng - vs)] (1) 

with Ecn = e 2 /[2Csjv(l — r 2 )]. The asymmetry pa- 
rameter b and the coupling parameter r < 1 are not inde- 
pendent, and r < min (6, 1/Vb). Eq. (Q~|i determines the 
charge stability diagram of the isolated NS system in the 
(vn,vs) plane. First, for a value vg imposing an integer 
number of pairs in S, say vg — 2, the charging number 
njv increases monotonously with vn- Next, consider a case 
where ng fluctuates, for instance vg = 1. For small r, 
as shown in Fig. 2a, tin is again a monotonous function 
of vn'- the sequence of charge states (nN,ng) as vn in- 
creases reads (0, 0), (0, 2), (1, 0), (1, 2), (2, 0), (2, 2), . . . (no- 
tice the oscillation of ng). The result is very different if r 
is large. In fig.|2jb), for vg = 1, tin increases with vn but 
in a non-monotonous way, the charge state sequence being 
(1, 0), (0, 2), (2, 0), (1, 2), (3, 0), (2, 2), etc. The correspond- 
ing charging staircases are plotted in the insets. 

One sees that the transition from (tin, 2) to (tin + 2, 0) at 
vn = tin + 1 "skips" the charge state njy + 1 in the grain. This 
signals a negative effective charging energy in N which over- 
comes the Coulomb repulsion. After increasing by two units, 
n n decreases by one unit, yielding a negative differential ca- 
pacitance (NDCA) Cdiff = C g N{driN / <1vn) at half-integer 
values of vn- Strikingly, the total number of steps, positive 
or negative, is doubled with respect to the usual case. Both 
charge skipping and NDCA occur above the dotted line indi- 
cated in Fig. 4 displaying a (b, r) diagram. From the above 
charging energy, an effective attractive potential U < can 
be estimated for v s = 1 as U = E c (0,2) + E c (2,0) - 
2£ c (l,2) = e 2 (l - 2rVb)/[CxN(l - r 2 )]. The necessary 
condition for the occurrence of a negative charging energy is 
thus 2r\fb > 1. 

To further analyze this possibility in an open NS system, let 
us consider the full Hamiltonian: 

H = E C + J2 £k C /U<x C kR,* + J2 £ 1 °In,* c qN,o 
kcr qa 

+ E Tfe >« - \ \ n s + 2 > (»s I + H.c] , (2) 

kqa 
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FIG. 2: Charge stability (or honeycomb) diagram for 6 = 1, r = 0.2 
(a) and r — 0.8 (b). The insets show the charging curves for N, 
taken along the dotted line. In case (b) the charging staircase (inset) 
exhibits charge skipping effects. 

where k (q) denotes electron states in the normal reservoir R 
(grain N), and the Coulomb interaction Ec is given by eq. (Q]). 

The total charge in N is expressed as tin = c^ N a c q N. a - 

go- 
Assuming constant densities of states in N and R, the single- 
electron transition rate from R to N is given by r( +1 ) = 

[5E ( c 1} /e 2 R N ] [exp (SE^ /keT)-!]- 1 within the golden 
rule approximation, where Rn is the tunnel resistance. 

Considering first the case of small Ej -C Ecs — 
e 2 /[2Cj;s(l — r 2 )], we perform a T-matrix calculation of 
the transition rates from (0, 2) to (2, 0) (close to vn = 1) 
and from (2,0) to (1, 2) (close to vn — 1.5). For the first 
transition, we take into account three possible configuration 
paths involving higher-energy states: (0,2) — ► (1,2) — ► 
(2,2) -> (2,0), (0,2) -> (1,2) - (1,0) (2,0), and 
(0,2) -> (0,0) -> (1,0) -» (2,0). For the second transition, 
only one excited state is involved: (2,0) — > (1,0) — > (1,2) 
and (2, 0) — > (2, 2) — > (1, 2). The shape of each step is cal- 
culated at finite temperature by solving the master equation 
governing the dynamics of the probabilities p(0, 2), p(2, 0) 
for the positive step and p(2, 0), p(l, 2) for the negative one. 
The master equation reads p(a) — T b ^ a p(b) — T a ^ b p(a) with 
p(b) = 1 — p(a) for the states a, b involved in the transition. 
Here the probabilities of other states are neglected, which is 
justified close to vn — 1 or vn — 1.5 and if the steps are 
sufficiently narrow. The resulting steps are shown in Fig. [3] 

As a result, a positive step +2e (where the charge number 
tin = 1 is skipped) and a consecutive negative step — e are 
stabilized. Notice that contrary to the usual case where all 
transitions between n and nil are real and obey the same 
master equation [19], here the rates are of higher order and 
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FIG. 3: Charging staircase : r = 0.8, b = 1, Ecs = 
E C n, k B T/E C N = 3 ■ 10" 2 , Rn/Rk = 10; (full line) 
Ej/Ecn = 0.5 ; (dotted line) Ej/Eon = 2. 




FIG. 4: Phase diagram in the b, r plane. The unphysical grey region 
is excluded. Charge skipping and NDCA occur above the dotted 
line, from bottom to top : Ej/Ecn = , and 1,2,4 (adiabatic 
calculation). All other parameters are the same as in fig. 3. 

the virtual states involved in one transition become real states 
(with first order rate) for the next one. A full treatment of all 
processes is beyond the scope of this Letter. 

Let us now turn to the case of a large Josephson energy, 
Ej > Eqs- The adiabatic assumption 01611 . setting the phase 
difference to across the JJ, allows to solve the Hamilto- 
nian (2), neglecting the normal electron tunneling term. The 
resulting adiabatic Hamiltonian H aa - = Eq — Ej cos <f> de- 
scribes a Cooper pair box with an effective gate voltage, which 
is an adiabatic function of njy. In the tight-binding limit 
Ej/Ecn 3> b, assuming that the junction dynamics is con- 
fined to the lowest Bloch band, one obtains the sum of the N 
dot charging energy and the adiabatic Bloch band energy : 

Had = E C n(1 ~ r 2 ){n N - v N ) 2 

r 

-A cos[7r(^s j=(n N - v N ))\ , (3) 

VO 

where the bandwidth is given by ffl 

A„ = l 6 ^E c „(^) 3 Vv^7^. ,4, 

The second term in H a d represents a nonlinear screening po- 
tential acting on the charge in N. The offset i/$ controls the 
phase of the cosine term, and an appropriate choice (for in- 
stance us « 1) achieves a negative curvature of H a d, viewed 
as an effective charging energy Eq£ for the gauged charge in 
N, un—^n- The required condition reads {it 2 /2)r 2 Ao/[6(l — 
r 2 )] > 1, yielding the lines in Fig. 4. Clearly, a large Ej 



puts a stronger constraint on the coupling capacitance Co, and 
requires larger values of r than for small Ej. The shape of 
the charge skipping and negative steps is then calculated, us- 
ing a master equation based on charge states = 0, 2 or 
iin = 2,1, respectively. The adiabatic transition rates are 
given by T ad = [SE e J^/e 2 R N ] [exp (SE e J^/k B T) - l]- 1 . 
The corresponding steps (Fig. 3) are less pronounced than in 
the small Ej case. 

Searching the optimum regime must account for the fact 
that for too large r values, the system behaves like one sin- 
gle island and its energy no longer depends on the location 
of the charge. The Coulomb blockade requires temperatures 
much smaller than the energy difference between two charge 
states. An optimal r is close to 0.75 (with 6=1) for small 
Josephson energies, and the requirement for Coulomb block- 
ade reads ksT < Ecn/^- The value r = 0.8 was used 
in the steps calculations. A temperature of T ~ 30 mK 
and a typical charging energy of Eqn ~ 0.1 meV were 
used, yielding C N = C j ~ 2 £F and C = 4C N = 8 fF. 
The tunnel resistance Rn is chosen as Rn/Rk — 10 with 
Rk = h/e 2 25.8 kil. This yields a bare single-electron 
tunneling rate of T ~ 10 9 s _1 . In presence of the S island, the 
effective tunneling rates are 10 7 s _1 (njv decreasing from 2 to 
1) and 5 • 10 3 s _1 (tijv increasing from to 2), respectively. 
With the above parameters, the maximum attraction \ U\ is of 
the order of 50/ieV. 

Notice that such a negative charging energy cannot render 
the dot superconducting, because it concerns the energy re- 
quired to add one or two electrons on the dot, rather than a 
true attractive potential felt by all electrons near the Fermi 
level. The effective "negative-U" potential manifests itself 
only when the dot is weakly coupled to a reservoir such that 
its charge can fluctuate. This is similar to a single "— U" cen- 
ter weakly hybridized with a normal bulk metal [9]. Here, pair 
fluctuations with the reservoir are assumed to be incoherent. 
On the contrary, a very small tunneling term Tjys between N 
and S would allow for a true phase coherence between states 
njv, tin +2, thus a proximity effect in a quite unusual regime, 
where T NS <\U\. 

Let us briefly discuss the issue of phase coherence in the 
Cooper pair box. Charge skipping only requires pair tunnel- 
ing between the superconducting reservoir and the S island in 
order to screen the repulsive interaction in the normal grain. 
No phase coherence is needed, as shown by our calculation 
in the small Ej case. Moreover, even in the large Ej case, 
charge fluctuations in N should strongly react back upon S and 
reduce the phase coherence. A full treatment goes beyond the 
adiabatic approximation [16]. One can anticipate that large 
fluctuations in the phase <f> renormalize Ej downwards, mak- 
ing the small-i?/ case generic. 

We now propose a scheme for detecting the non- 
monotonous charging of the N grain. SETs, or point contacts 
1I20I1 provide very sensitive detection of the local change in the 
electrostatic potential. In double-dot setups with weak mutual 
coupling, the potential variations in each dot can be measured 
by a different neighboring point contact l2ltl . In the present 
case, placing a point contact close to N does not measure Sn^, 
but instead SVn — (C _1 )7vw(e^«Ar) + (C^ 1 )Ns(eSns) = 
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e[5n N + rVbSn s ]/[C^ N (l - r 2 )]. If 2rVb > 1, dou- 
bling of the number of steps can be detected, but not the non- 
monotonous charging curve. To access the latter, it is suitable 
to measure SVs — e[ry/bSnN + bSns]/[C^N(^ — r 2 )] as 
well, with a second point contact close to S, and reconstruct 
SriN = Csn [5Vn — rSVs / Vb]/e. The parameters Ceat, r, b 
can easily be measured from the stability diagram obtained 
in the normal state in the presence of a very weak tunneling 
between N and S lfl8ll . Notice that the tunneling rates calcu- 
lated above are much reduced compared to the bare single- 
electron rate T. Therefore the use of point contacts permits, 
non only a time-averaged [21], but even a time-resolved and 
directional [22] detection of the charge variations in N and S. 
On the other hand, cross-correlation shot noise measurements, 
as in ref. 12311 . would require higher currents. In practice, a 
possible setup inspired by ref. 112 ill is proposed in Fig. Q] It 



involves a superconducting strip with a Cooper pair box, cou- 
pled laterally to an InGaAs/AlGaAs 2DEG, suitably tuning 
the Schottky barrier present at the interface between the su- 
perconductor and the 2DEG. Notice that the geometry of Fig. 
1 could be modified, including drain and source such as to 
allow transport through N |0] • 

In conclusion, we have proposed a mechanism inducing a 
controllable negative charging energy, thus attractive correla- 
tions, in one or several metallic dots. We believe that such an 
effect would be useful in view of more complex nanoelectron- 
ics devices. The authors are grateful to T. Martin, M. Fogel- 
strom, and G. Johansson for useful discussions, and S. Ander- 
gassen for careful reading of the manuscript. D. F. and A. Z. 
were partially supported by the contract AC NANO NR01 14, 
and C. H. was supported by the Swedish Research Council 
(VR) under grant 621-2006-3072. 
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